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ABSTRACT 


This fchesis presents a study of state feedback model reference euiaptive control 
technique using Lyapunov's second method for the design of power system stabilizer (PSS). 
The reference model incorporates the desired chareu:teri8tic8 of system performance. 

A numerical power system problem has been considered and an adaptive PSS 
designed for it Detailed simulation study of the behaviour of the controlled system has 
been done at various operating conditions and system strengths to establish the efficacy of 
the procedure. 

Since all the state variables are not available in practice, €uad since even the system 
parameters are not known, reconstruction of the state variables by the usual procedure is 
not possible. Here, smGe)tthe adaptive system works properly, the controlled system tends 
towards the reference model due to the adaptation process, it was felt that observers based 
on reference model parameters and plant outputs may work satisfactorily. Detailed 
numerical experimentation has established that a third order observer, which uses only the 
angular velocity signal &om the generator, cascaded with the state feedback adaptive 
controller gives highly satisfactory results. 
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CHAPTER 1 
INTRODUCTION 


To improve system trcmsient stability, fast acting high gain AVRs with static 
excitation systems are used with modern alternators. But, they have deleterious effect on 
dynamic stability, resulting, often, in lightly damped and even sustained 
oscillations, typically in the range of 0.2 to 2 Hz. Thus it can be assumed that the modern 
voltage regulator introduces negative damping. To offset this effect and to improve system 
damping, artificial means of producing damping torques have been developed. These means 
essentially produce torque components in phase with generator rotor speed and involve use 
of an auxiliary circuit in the excitation control loop. The output of the auxiliary circuit ia 
added to the main voltage feedback of AVH. These circuits have been traditionally called 
as power system stabilisers (PSS). Their input is a signal from rotor speed, accelerating 
power or generator bus frequency. 

Various methods have been reported in the literature for the design of PSS. Earlier 
designs were baaed on classical control theory and the recent works are based on moderil 
control theory. The following paragraphs give a brief review of the important works done in 
PSS design. 


1.1 Adaptive PSS 

tn tMe Asign method, the PSS structure and parameters, once designed, will not 
change sdth timror operating conditions or system configuration. Such PSS will give good 
results if tW curating condition does not change much from that used for PSS design, 
An extentfi^ |^et of flAthods falling into this category is available in the literature; see for 
example [22,^, pnljf »few of them are discussed, here. 

Extensive study in this field started with the paper <f dc Hello and Concardia[l] 
which explained the problem of dynamic stability in terms of the concepts of synchronous 
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and damping torques and proposed design of PSS using classical frequency domain 
techniques. Bollinger et. al. [2] used the root locus technique to shift the eigen values, 
corresponding to rotor oscillations, to suitable locations Hamdan and Hughes [3,4] used 
multivariable frequency response plots for design and Euialysis of PSS. Based on linear 
regulator theory, Yu et. al. [5] have designed the stabiliser as an optimal state feedback 
controller, which mimmises a quadratic performance index (PI). The resulting controller, 
in this method, improves the closed loop performance, but the inaccessibility of some of the 
states for feedback creates a problem in practical applications. Suboptimal controllers with 
the system represented by lower order models, with only measurable states as state 
variables, have been developed by Yu and Siggera [6]. 

Pole assignment techniques have been considered by many authors for the problem 
of power system stabilization. This approach provides a design method lor PSS with the 
control objective of assigning the dominant closed loop eigen values to specified stable 
locations using either state feedback or output feedback. Prom a practical view point, pole 
assignment using output feedback is more important. Padiau et. al. [7] have presented three 
methods for the design of PSS by pole assignment with output feedback. Srinivas and 
Rameir [8] have designed first and second order PSS using pole placement with output 
feedback. Chow and Sanchez-Oasca [9] have presented four methods for PSS design, 
namely, full state feedback, output feedback by projective control method, full order 
observer and dynamic compensator by the projective control method with main emphasis 

on frequency response charecteristics of PSS transfer functions. 

' I . 

A ijMlstice, a power system will be very large consisting of hundreds of generating 
units, tts Buch a large system will be very difficult, if it is represented in terms 

''' tpi 1 

of the svuvalent circuits of each and every machine. Undrill et. al, [24] have 

proposed Ibif the analysis of large power systems by separating the whole power 

system intO'^^^parts, namely, the study system and the external system. The study 
systSih consist of & few machines of extreme importance and it is represented exactly. The 
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external system consists of the remaining machines, which are very large in number but of 
secondary concern These machines are represented by their dynamic equivalents 

Rogers and Kundur, in the discussion of [25], recommended the dynamic reduction 
technique of Gremond and Padmore [26] for the analysis of large power systems. In this 
technique, each coherent group of machines is represented by an equivalent single machine. 
This technique is very effective in reducing system size. 

Wong et. al [28] have developed an algorithm for the analysis of dynamic stability 
of very large power systems consisting of more than thoussmd machines. This method is 
very effective in determining the interarea stability. Here the exciter modes are masked 
and a search is made for electromechanical modes. This algorithm is based on the 
essentially spontaneous oscillations in power systems, developed by Byerly et. ed. [27]. 
Unlike the early studies, where D Q (direct axis and quadrature axis) variables are used, 
here, P N (positive sequence and negative sequence) variables are used. 

Uchida and Nagao [29] have proposed an S matrix method for the determination of 
eigen values of large electric power systems. In this method, the A matrix of the plant is 
transformed into the S matrix and the dominant eigen values are identified. This method is 
very fast due to the spcu-sity and structural uniformity of S matrix. 

1.2 Adaptive PSS 

A stabiliser designed using conventional control theory at one operating point may 
not give satisfactory performance at other operating points. In such situations, adaptive 
PSS are to be employed, which automatically sense and correct themselves whenever 
significant changes in system parameters occur. All the major adaptive control techniques, 
namely, model reference ade^tive control (MRAC), seif tuning regulators (STR), gain 
scheduling schemes £Uid variable structure control have been considered for adaptive PSS 
design. A brief review of the major works in thm areas is given below. 
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L2.1 Self Tuning Regulator (STR) Techniques 

Design of adaptive PSS as a STR has become feasible because of developments in 
microprocessor technology. These are digital adaptive control schemes and involve 
par£nneter estimation and control. In explicit STR, parameter values are explicitly 
obtained and used for generating control signals. In implicit schemes, parameter estimation 
is implied but is not done explicitly. Recursive least squares procedures are the most 
popular for identification 

Ghoeh et. al. [10] summarise the features of MRAC and four versions of self tuning 
controllers (STG) for ade^tive PSS. In [11], the same authors have designed an adaptive 
PSS for large disturbances using recursive least squares (RLS) technique for parameter 
identification. By testing the PSS on an experimental micromachine system, it has been 
shown that this self tuning PSS gives better results than those with a fixed parameter 
stabilizer. Kanniah et al [12, 13] have used an implicit type STO, with RLS technique for 
parameter identification. Simulation results in [12] show that the dual rate STC, where two 
different sampling rates, one for identification and the other for control, are used, performs 
better than the single rate STC. In [13], the hardware and software used to implement dual 
rate STO for controlling a laboratory alternator system, are described. Xia and Heydt [14] 
have developed an STG approach for PSS design for a generator connected to an infinite 
bus through a long transmission line. The salient points of the controller are its adaptive 
nature and application of a least squares formulation for parameter estimation. The power 
system is stable over a wide range of operating conditions. Ohandra et. al. [15] have 
implemented a self tuning excitation controller performing the functons of AVR and PSS 
along with an interactive man— machine interface using two microcomputers. 

1.2.2 Gain Scheduling Schemes 

In this scheme, a range of operating conditions, which are encountered in practice, 
are taken into consideration. Extensive ofiOine design and simulation effort are needed to 
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obtain a suitable gain scheduling scheme. So, this scheme is very time consuming. 

A simple ad€ 4 }tive PSS using gain scheduling approach has been developed by 
Brown Boveri company [16]. For discrete values of operating points and system strengths, 
suboptimal gains are calculated apriori and stored in the memory of a microcomputer in 
the form of a three dimensional lookup table. Depending upon the operating point, they are 
retrieved for PSS tuning. Sen Gupta [22] used the properties of tensor invariance to 
isolate sources of positive damping and synchromzing torques. Synchronizing and damping 
torques are exactly quantified. The PSS developed cancels the negative damping generated 
by AYR, and aug: ments the positive damping part of its contribution. Ashok Pal [17] has 
designed PSS for large disturbances. An extended Kalman filter attached to a detection 
scheme is used for identification. Using the results of identification and after detailed 
digital simulation, two lookup tables, containing optimum PSS gains for different operating 
points, are formed. Using these lookup tables and using linear interpolation, PSS 
parameters are obtained for any given operating point. 

Bandyopadhyay and Prabhu [30] have proposed a new approach for the design of 
adaptive PSS using gain scheduling technique. Five different constant gain PSS are 
designed for five operating conditions. All these five PSS operate in parallel. At any 
operating point, the control signal for stabilization is obtained as a weighted sum of the 
outputs of the individual PSS. The weighting factors are determined using a simple 
algorithm. This adaptive composite PSS (AOPSS) performs satisfactorily over a wide range 
of operating conditions. 

Recently, Madhu and Prabhu [31] have proposed a method for the design of robust 
adaptive composite PSS (RAOPSS) making implicit use of gain scheduling schemes. In this 
method, a set of PSS is designed. Depending upon the operating condition an appropriate 
PSS is selected from the set using a simple algorithm. It is shown here that a maximum 
three individual PSS are sufficient for satisfactory performance of the system. The 
algorithm used here is simpler and easier than that used for the design of AOPSS. 
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1.2.3 Variable Structure Schemes (VSS) 

VSS essentially involve nonlinear state feedback control. The control law involves 
switching, whenever the system state crosses a 'awitching hyperplane' in the state space. 
Since the structure of the feedback changes depending on the position of state point, these 
schemes are called variable structure schemes. Once a state point reaches a switching 
hyperplane, it is made to travel on the hyperplane, during control Thus, as soon as a 
disturbance changes a state point from its equilibrium vedue in a system with regulating 
control, the controller brings it back to the equilibrium point by making it move on the 
switching hyperplane. Motion of the system on the switching hyperplane is called 'sliding 
mode of operation An advantage of these controllers is that control action takes place 
even in the presence of parameter variations. In many cases, the physical realisation of such 
a controller can turn out to be very simple. Switching action permits the use of limiting 
values of control input magnitudes thus leading to high speed of response. The only 
drawback of these schemes is the excessive 'chattering' of the controller which is 
undesirable, in practice, &om the point of view of both controller and plant hardware. Hsu 
and Chan [18] have attempted PSS design as a variable structure stabilizer. Serious work in 
VSS for power system stabilization has not been done. Strictly speaking, VSS should be 
looked at as adfq}tive controller, since it does not attempt to change the controller on the 
basis of variation in plant parameters. However, it is possible to develop VSS in the MHAO 
setting as is done by Ambrosino et. ai. [32]. 

1.2.4 Model Reference Adaptive Control (MRAC) Technique 

Irving et. al. [19] described MHAO for generator voltage regulation based on 
maintaining an unconditionally stable adaptive loop using Popov's hyperstability theory. 

Ibrahim and Kamel [20] have proposed a model reference adaptive PSS (MRAPSS) . 
They claim that this PSS improves the dynamic stability of power system by effectively 
increasing <ihe damping torque of the synchronous generator in the system. They have used 
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Lyapunov technique to assure stability of the system with PSS. In this work, they have 
chosen a second order transfer function for the power system, with torque as output and 
speed as input, with no zeros. But, the transfer function is in fact a third order one with 
three poles and two zeros. They have tried to design a first order PSS. But, a first order 
PSS is not at all acceptable for practical purposes. They have compared a component of 
electrical torque, which is proportional to the q— axis component of voltage behind transient 
reactance, of the plant with that of the model. But, this component of torque is not at all 
physically available separately from the total electrical torque. Simulation results are 
wrong, since even the calculation of well-known coefficients k^, k^, k^, k^, k^ and k^ is 
wrong. So, this scheme, under no circumstances, can be considered for the design of PSS. 
Qhandakly and Idowu [21] have presented a decentralized MRAC scheme for the design of 
PSS and a means for coordinating the generator unit excitation and governor control loops. 
System stability is assured by deriving the adaptive control law from Lyc^unov energy 
function. 

1.3 Scope of the Present Work 

The present work is concerned with computer aided design of MRAPSS, which 
performs satisfactorily over a wide range of operating conditions. System stability is 
assured by deriving the adaptive controller from application of Lyapunovas second method. 
In this scheme, the state variables of the power system under consideration (the plant) are 
to follow those of an explicitly specified reference model, which is designed to have 
desirable performance charecteristics. In general not all the state vaunables are available for 
feedback. So, the unavailable states are recozistructed using reduced order observers and 
the effectiveness of MEAPSS is investigated. 

L4 Outline of the Rest of the Thesis 

Cbspter 2 describes the power system considered. A brief discussion on modelling of 
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the power system for dynamic stability studies is also ^ven It is followed by a brief survey 
of the PSS input signals used in the traditional PSS, ia the literature. Then, the FSS input 
signals used in the present work, are ^ven. 

In chapter 3, a brief introduction to model reference adaptive systems (MBAS) and 
their design methodologieB are given. The PSS design problem is posed as a MBAS design 
problem. Lyapunov's synthesis technique is used for designing the adaptive controller, 
which ensures system stability. A numerical problem is considered. First, assuming all the 
states to be available, a model reference adaptive controller with parallel structure and 
having full state feedback is determined and the resulting system is simulated over a wide 
range of operating conditions. Since all ^^states are not available, in practice, reduced 
order observers are next designed for reconstructing the unavailable states. Design of the 
observers needs knowledge of the plamt parameter matrices. These matrices, however, are 
not known. We can use €ui adaptive observer along with a state feedback adaptive 
controller. However, what happens when the two systems, namely, the adaptive observer 
and the adaptive controller, designed seperately, are used together, is not known. 
Furtb^nnore, the resulting system complexity may discourage its use in practice. So, 
instead of the above approach, we have used the parameter matrices of the reference model 
for state observation. Since, in steady state, the plant with the feedback controller is 
expected to converge to the reference model, it is hoped that, asymptotically the state 
observation leads to the actual state values of the plant and the adaptive system would be 
stable, just as the full state feedback MBAOS designed using Lyapunov's direct method. 

Assuming that one, two and three state variables are unavailable at a time, first, 
second and third order observers are respectively constructed, and the system responses are 
examined for all the three cases over a wide range of operating conditions. Good results 
have been obteuned suggesting that the ^proach taken up is effiicacious. 

Chapter 4 contains conclusions and suggestions for further work, 


CHAPTER 2 

MODEL OF POWER SYSTEM 


High gain, fast acting modern voltage regulators often lead to sustained low 
frequency oscillations in power systems. In order to damp out these oscillations and assure 
dynamic stability, an auxiliary controller is used in the exciter control loop 

In this chapter, the model of the power system considered for the study of dynamic 
stability is presented. 

In its simplest form, a power system consists of a synchronous generator connected 
to an infinite bus through a transmission line, In the literature, various state space models 
of this system having varying degrees of complexity and accuracy are avulable. In the 
classical model, the synchronous machine is represented by a voltage source behind the 
machine direct axis transient reactance. Though this model is very simple, it does not 
account for the demagnetising effect of armature reaction. 

The papular Hef&on—Philips model has been widely used by many researchers for 
the study of dynamic stability. This model is almost as simple as the classical model and 
offers a good physical understanding of the synchronous generator while taking the 
demagnetising effect of armature reaction into consideration. It neglects the fast dynamics 
associated with the damper windings since these are not of importance in the usual 
dynamic stability studies. Heffron— Philips model and the slightly more complex model* 
[33] have, perhaps, become the standard models for dynamic stability studies. There are, 
however, problems for example those associated with subsynchronous resonance (SSR), 
where more detailed models would be needed. The models used for dynamic stability 
studies are all linear models defined etfound the system operating point. 

In the present work, we have confined our attention to the He&onr-Philips model. 
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2 A Heffron-Philips Model (State Space Representation) 

In the development of the model, the following assumptions are made: 

1. The effects of damper windings are neglected. 

2. Stator winding resistance is neglected. 

3. The speed voltage terms are dominant over the rate of change of flux linkages 
terms. 

4 Balanced conditions are assumed and saturation effects are neglected. 

5. The electrical network is considered to be in quasi— steady state. 

The power system considered is shown in Fig 2.1. It consists of a single machine 
with an AYR connected to an infinite bus throu^ a transmission line. The Heffron— Philips 
model of the system around an operating point is represented in block diagram form in Fig 

22 ,&^ 3 ' 

The Heffron-Philips model of the power system can be expressed in state space 
format as 

x=Ax+Bu (2.1) 

y=Ox 

where 


< 

II 

X 

"A 



C2..3: 


■ 0 

1 

0 

0 ■ 



A — 


0 

“2 

0 



A— 

0 

0 

“3 

“4 





0 


"V ■ 



B=[ 

0 

0 

0 

b]‘ 


C2-’5j 


The output variable y is defined to be that which spears as PSS input. Many 
different PSS input signab are used in practice. The matrix 0 thus depends on the input 



I I 3'do| 

- I V.A° ^5W;A , 


Fig 2.2 Block Diagram Of Power Bysfcem 
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signal used by the PSS in a given situEition. 

A brief derivation of the model is given in Appendix 1 along with definitions of 
coefficients ^ ’ ^5 ' *6 ' ^ 

2.2 Input Signals to PSS 

In PSS design, proper selection of input signals is important. Feedback of rotor 
speed deviation appears to be proper, since the objective of PSS is to provide adequate 
damping of electromechanical oscillations. Speed input signal has been used by many 
researchers. Other signals, such as terminal voltage, armature current, active and reactive 
powers, accelerating power and frequency can also be used as input signals to PSS 

Larsen and Swann [23] have analysed the effectiveness of speed, power and 
frequency signals for PSS> They have used frequency domain arguments based on classical 
control theory. 

Anwar [35] has examined the suitability of PSS input signal in terms of stability 
robustness of the system, i.e., in terms of regions in the P— Q plane, where the system 
remains stable. In this study, he has used a state space description of power system. Here P 
is the active power and Q is the reactive power of the generator From the practical 
experience at Ontario Hydro, Canada, Kundur [36] has given a brief discussion on the 
efrectiveness of various input signals. Madhu [37], utilising the idea of robustness, clearly 
showed the superiority of two input PSS, with speed and acceleration as inputs, over single 
input PSS. 

In the present work, we use the full state vector as PSS input signal. In practice, all 
the state variables considered are not directly aveulable for measurement and will have to 
be obtained throu^ some kind of observation scheme, These matters will be taken up in 
the sequel. 



CHAPTERS 

MRAC TECHNIQUE FOR PSS DESIGN 


In the first chapter, the importance of PSS in damping out low frequency 
electromechanical oscillations in power systems has been described. It has also been 
indicated that a fixed parameter stabilizer designed at a particular operating point may not 
give satisfactory performance when the operating point changes. In practice, the operating 
point, which depends upon load conditions, changes continually. Therefore, should 
adapt itself to the changing load conditions. Unlike fixed gain or fixed parameter 
stabilizers, the adaptive stabilizer determines a new set of control parameters as changes 
occur in system configuration and load levels, thereby ensuring that the controller 
parameters are effice^ious for a wide range of operating conditions. The model reference 
approach for PSS design has the definite advantage of stability assurance which is lacking 
in self timed regulator design. Furthermore, higher adaptation speeds can be obtained with 
MBAC, [38]. 

A design technique, which involves the application of MRAC concept for the 
purpose of designing a PSS placed in the excitation control loop is described in this 
chapter. In this scheme, a reference model into which the design specifications are 
incorporated is introduced and the states of the plant are forced to follow the states of the 
specified reference model. Adaptation signals are derived for the exciter input to assure 
overall system stability as well os the desired system performance. This design scheme is 
verified on an example synchronous generator connected to an infinite bus through a 
transmission line. Since all the states of the generator are not always available, numerical 
experimentation is done to know what happens if the unavailable states are reconstructed 
using state observers. By simulating the overall system over a wide range of operating 
conditions and for different system strengths, it is shown that a third order observer which 
reconstructs states and using the state gives satisfactory design of PSS. 


14 


3.1 Introduction to MRAS 

MRAC systems are important since they lead to systems with a fairly high speed of 
adaptation, as there is no need for estimation of plant paremietem. In MBAC, the aim is to 
m€Jce the states of the unknown plant approach asymptotically the output of a given 
reference model, which is part of the control system. The di^erence between the states of 
the reference model and those of the adjustable system is used by the adaptation 
mechanism either to modify the parameters of the adjustable system (called parameter 
adapiaitori) or to generate an auxihary input signal (called signal synthesis adaptation) in 
order to minimise the difference between the states of the adjustable system and those of 
the model. The basic scheme of MRAC is given in Fig 3.1 

In terms of controller structure, MRAS are classified into three types; 

1. Series MRAS: 

Here the adjustable system and the reference model are in series, Fig 3.2. 

2. Series Parallel MRAS: 

There are two cases here. 

i. The reference model is in two parts, one in series with the adjustable system and 
the other in parallel, Fig 3.3 1 

ii. The adjustable system is in two parts, one in series with the reference model and 
the other in parallel, Fig 3.3.2. 

3. Parallel MRAS: 

This is the most commonly used scheme. Here, the reference model is in parsdlel 

with the adjustable system, Fig 3.4. 

We have confined our attention to the parallel MRAC structure in the present 


work. 
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Fig 3.2 Series MBAS 
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Fig 3.3.1 Series Parallel MBAS 


Fig 3.3 2 Series Parallel MBAS 



Pig 3.4 Parallel MBAS 
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3.2 Design Methodologies 

Most of the techniques proposed for design of MRAG reduce the adaptation problem 
to that of determining the controller structure and finding how the controller pareuneters 
are to be adjusted, that is , obtaining a parameter adaptation law to guarantee global 
asymptotic stability. There are three important techniques for the design of continuous 
time model reference adaptive controller, namely, the MIT rule and related methods, the 
Lyapunov approach and Popovas hyperstability i^proach. 

3.2.1 MIT Rule 

The MIT rule does not assure system stability. A good compromise is made between 
stability and speed of adaptation. 

Here, the input r(t) of the system is e^pUed to appropriately designed reference 
model to generate the desired output. The integral of the square of the error between this 
desired output and the output of the plant is minimised. 

The performance index (PI) is given by 

J = (3-1) 

0 

where 

y is the output of the reference model and 

y is the output of the plant (adjustable system) 

It is assumed that the speed of adaptation is low and the system parameter values 
are always in a sufficiently small neighbourhood of model pareuneter values. 

Since the minimisation is to be performed online, the performance can be evaluated 
over a sufficiently small time interval At. Over At, PI is 

= 3 II ll' 

Here b is the adjustable parameter. 

J^(b) is dependent on t and is to be minimised with respect to b by using a search 
techni^e, after it has attained a steady state value. The search techniques are 
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implemented using sensitivity functions, which must be generated online 
One of the search techniques to generate the adaptation algorithm is the gradient 
method. The standard gradient search algorithm may be used to generate estimates of the 
parameters. 

b(t) = -B ^ = B ^ o(t) (3.3) 

db db 

where 


B is an arbitrary positive constant called the adaptive gain, 
e = (yjt)-y(t)] (3.4) 

and ^ is the sensitivity matrix. 

The equation (3.3) gives the MIT rule. 

The MIT rule was popular because of its simplicity in practical implementation. 
But this must be used with caution, inview of possible instability. Furthermore, some other 
drawbacks are; 

1. The initial difference between the parameters of the reference model and those of the 
adjustable system is assumed to be small. 

2. The adaptation speed is also assumed to be less. 

3. There is no indication given regarding the selection of the adaptation gains in order to 
assume the convergence of the adaptation process and the stability of the MBAC B3rstem. 

3.2.2 Lyapunov Approach 083 

In order to achieve good design, it is to be made sure that the system along with the 
controller is stable. Stability can be ensured by using Ly^unov approach or hyperstability 
approach in designing controllers. 

In Lyapunov approach, the first step is to determine the differential equation which 
describes the error between the states of the model and that of the plant. Then the 
puameter adjustment equations are obtained as the conditions to assure the error 
differential equation, stable. ,ln order to do this, a positive definite Lyapunov function isi 
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chosen. Then, the ad£ 4 )tive mechanism equations are derived so as to cause the time 
derivative of the Lyapunov function negative aemide£nite This negative semidefimteness 
of the Lyapunov function ensures system stability and the error is also made to go to zero, 
as t -» OD. 

The main difficulty with the Lyapunov ^proach is that the entire state vector must 
be available for measurement and this is not always possible. Another drawback is that the 
Lyapunov design rule may not be applicable to cases where the plant parameters can not 
be directly adjusted. In cases, where a Lyapunov function can not be found, the 
hyperstability approach can be used. 

3.2.3 Hyperstability Approach 

Popov [39] introduced the concept of hyperstability as a natuisJ extension of 
absolute stability. Landau [40] proposed application of this concept for the design of 
adaptive control systems. 

For application of hyperatability conditions, the MRAG system (MEAOS) is 
expressed as a standard system (Fig 3 5), with a linear time invariant operator in the 
forward path and a passive non linear operator in the feedback path. The non linear 
operator satisfies Popov Integral Inequality, [38] The system is then asymptotically 
hyperstable if the transfer function matrix of the linear block is strictly positive real, [38]. 
To apply this approach to MRAOS, it is necessary to express the error equation in the 
standard form. Hyperstability conditions enable determination of feedback block, N 
(Fig.3.5) which constitutes the adaptation law. 


n -> 

L 

V 

? 



u 

N 

^ 


Fig 3.5 Standard Feedbcu^ System 
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Landau’s Procedure for the Determination of MRACS, [38]: 

Step 1 ; The MRACS is transformed into an equivedent standard feedback system 

composed of two blocks as shown in Fig 3 5, which is always possible. 

Step 2.: The solutions for the part of the adaptation laws which appear in the feedback of 

the equivalent system are found such that the Popov integral inequality, 

T 

/ u‘ V dr > (3-6) 

is satisHed. Here, u and v are as shown in Fig 3.6 

Step 3 : The solutions for the remaining part of the adaptation law, which appears in the 
feedforward path, are found such that the transfer function matrix of the feedforward path 
is strictly positive real, This ensures global asymptotic stability of the entire system. 

Step 4 : Return to the original MRAOS in order to specify the adaptation law explicitly, 
that is, the structure of the mechanism. 

Since this approach requires the use of derivatives of error, when only some of the 
states of the model and adjustable system are available, this gives a significant diHLculty. 
The Lyapimov approach and the hyperstability are equally efEcient, but in many cases, 
hyperstability approach may be easier to apply than Lyapunov approach. 

In the present work, we have used state feedbevck MRAO method with Lyapunov 
synthesis technique for the design of PSS. 

3.3 Comparison of MRAC with STR 

MRAOS are easy to implement with a fairly high speed of adaptation. They can be 
used in variety of situations. MRAO can be used with continuous as well as discrete 
systems, whereas STE can be used only with discrete systems. MRAOS are mostly ^plied 
in deterministic systems, whereas STR are essentiidly stochastic systems. Recent research 
showed the equivalence of MRAO and STR methods. Well designed STR have the 
advantage of being able to adeq}t to any disturbances, MRAOS have the definite advantage 
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that the design ensures over^l system stability. This assurance is lacking in STR 


3.4 State Feedback MRAC Using Lyapunov Synthesis Technique 

The basic configuration of state feedback MRACS is shown in Fig 3 6 


-^I'ReFerence 
Model 


u 


Plant 


Adaptive 

Controller 


Fig. 3>6 Basic Configuration of State Feedback MRACS 


The state equation of stable reference model is 


X = A X + B r 
m m m m 


and that of the linear time invariant plant is 

X = Ax + Bu 


(3.6) 


(3.7) 


Here, x and x are n— dimensional state vectors and r and u are ni-<limensional 
' m 

control vectors. 

The state error vector is defined as 


e = x^-x (3.8) 

The differential equation for the error is ^ven by 

4 = A e + f (3.9) 

SI 

where 

f = (A^ - A)x + B^r — Bu , 

In order to drive the state error vector to zero by feedback synthesis using 
Lyapunov's direct method, a positive definite function of the form 




is selected as a candidate Lyapunov function where P is a symmetric positive definite 
matrix and h(^, is a functional expression for the misallignment between the reference 
model parameters euid the plant parameters. Misallignment parameter vectors (p, ^md ^ are 
defined in terms of the elements of the matrices (A^ — A) and (B^ — B) respectively. 

The control objective for a stable system is to specify a control law such that the 
first derivative of the candidate Lyapunov function, V is negative semi definite. This 
implies the stability of error system. If an adc^tation law is so chosen that V s 0 implies 
e(t) = 0, then by a theorem of La Salle, e(t) -* 0 as t oo. This means the controlled system 
is asymptotically stable 

Now let us have a state feedback of Kx and a feedforward of Su for the plant. In 
that case, the plant state equation, with the above control incorporated, becomes 

x=(A + BaC)x+BSu (3.12) 

For exact model matching, we need existence of matrices K andS such that 


BE =B^ 

m 

A + B„K* = 

m m 


(3.13) 


It is assiuned that both B and B are of full rank so that E is non-singular. The 


m 


equation, satisfied by the state error vector, e is 

e = A^e + (A^ - A - BEK)x + (B^ -BE)u 
It can be shown that the above equation can be expressed as 


(3.14) 

(3.14) 


where 


^ = [K — K(t)] and 


(3.16) 


For the error equation, we choose 


-1 


-1 
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where 

= rj > 0 and 

= r‘ > 0 (3.17) 

It can be derived, [41], that, if we have 

^ = "r.B ' Pex* and 
1 m 

^ = -rjB^Pe (u + Kx)*S‘ (3. 18) 

The stability of equation (3.14) in the (e, tp, space is guaranteed In the above equation, 
P is a positive defmite matrix satisfying 

= - Q, Q = Q‘ > 0 (3.19) 

The adaptive control law in terms of k(t) and S(t) is 

li((t) = r^B^Pex* 

E(t) = SPjB^Pe (u + Kx)‘S^ (3.20) 

It should be noted that (see [41]) stability is not assured in general over the whole 
(e, k, S) space. A necessary condition for stability is that the initial values e(tp), t^t^) and 
I!(tp) should satisfy 

v(to) = I [e*(to) Pe(tQ) + tr { (k* - K(t,))* r^(k* - K(t„) + (e’"‘ - ir‘(t,))* 

—1 —1 * —1 
(E* -ir\))}] < ltr[(E* )Vp* ] (3.21) 

Thus, it is necessary to start with sufficiently small values for ^(tQ), k *- k(tQ) and 

s^-s(V- 

3.5 Posing of PSS Design Problem as a MRAC Design Problem 

The state space model of the power system for the purpose of PSS design has the 

form 

X =: Ax + Bu (3.7) 

where A is 4x4 matrix, B is 4x1 matrix, x is 4x1 vector and u is the scalar input. 

The reference model is represented as 
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*m= Vin+ V (3-6) 

The state variables in the above equation are 

X = [Xj Xj X3 xj‘ = (3.22) 

\ \ *zn/ = C^mA "mA ^^mA W^* (3.23) 

The control matrices B and B are 

m 

B = = [0 0 0 b]‘ (3.24) 

The matrix A is given in equation (2 4). 

From the derivation of the power system model given in Appendix 1, it is clear that 
as the operating conditions defined by the complex power supplied by the generator 
changes, only the matrix A chuiges, the matrix B is unaffected. Thus, it is not necessary to 
have the feed forward term Su in the design as given earlier and the adaptive control 
becomes 

u = l(S:t)x + r (3.25) 

where u represents the control input to the plant. The closed loop system becomes 

X = (A + BlOx + Br (3.26) 

The matrix k is given by 

li(t) = r (B^Pe) X* (3.27) 

where 

r is a symmetric positive definite matrix and 

P is symmetric positive definite and satisfies the Lye^imov equation (3.19). 

A^ is to be chosen such that it ensures the asymptotic stability of the reference 

model. 

Selection of reference model - 

The matrix A of the reference model is selected such that the poles of the system 
m 

are placed in the permissible sector shown in Fig. 3.7. The power system chosen by ua is 
essentially the same as that of Madhu [37] and we have chosen the same sector as was 
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chosen by him for PSS design by eigen value placement. 


Fig. 3 7 Permissible Sector of Eigen Values 

3.6 Numerical Problem (Pull Stats Feedback) 

For the data given in Appendix 2, at rated load and UPF 



0 10 0 

-86.662 0 -128 72 0 

0 0 -0.3914 0.2 

675.3 0 -3142.2 -100 

The open loop poles of the system are 0.66 ± j9 63, —8.29 and —93.2113. 

We specify that should have eigen values at — 1 ± j8, —3.0 and -93.2113. 

We will determine A^ by determining the 1x3 matrix Q such that A + BG has the 
required eigen values. Then A + BG gives us A^. It is found that 

0 10 0 

-86.652 0 -128.72 0 

0 0 -0.3914 0.2 

-796.8 -62.1 -2066.6 -97.82 

where G = [-0.274 -1.242e-2 0.216 4.36e-4] 

As mentioned earlier, B is chosen to be the same as B. 


(3.29) 

(8.30) 
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Let be represented as 


A = 

m 




0 0 


m. 


m. 


m. 


a 


m. 


m. 


Let r be chosen to be diagonal 


m. 




(3.31) 


r = 


7i 0 0 0 

0 72 0 0 

0 0 Tg 0 


(3,32) 

0 0 o' 74 J 

Q is chosen to be diagonal, and equal to ql, where q > 0 is a real scalar. For the 
chosen Q, the Lyapunov equation (3.19) is solved for P. Let the last row of P be [p^ p^ 
pj. Then, we study the resulting design, after incorporating the adaptive feedback law 
given by equeU)ions (3 26) and (3.27) by extensive simulation, for various operating 
conditions, using the software package SIMNON. The equations for simulation are the 
plant equations, the reference model equations and the nonlinear feedback controller 
equations. These equations are given below. 

Plant : 


Reference Model ; 




Xj = a^x^ + a^Xg 

^3 = ^^+ V4 

*4 = Vl + Vs + V4 


(3.33) 


X = X 

“l “s 
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X = a X + a X 
“2 “l “l “2 “3 


x_ = a X + a X 
“3 ®3 ®3 “4 “4 


*«»+*« * ^ 3C + a X + br 

”^4 ”5 “1 “6 “2 + “7 “3 “8 ”^4 


(3.34) 


Feedback Qaina ; 


S = 

'‘4 = '^m'!'4*4 (®-36) 

= -b [ Pi(*i - V) + Pa(*2 - *m ) + P3(*3 - *m ) + P 4 W " V) ] 

1 / O 4 

(3.36) 


Control Input : 

u = kjXj + kjXg + kgXg + k^x^ + r (3.37) 

By experimentation, it was found that P = I and q = 10~^ resulted in adeq^tive 
controller design which not only gave stable response over a large range of operating 
conditions, but also exhibited good damping in the response. Some representative plots are 
shown in Figs. 3.8 to 3 13. Here the operating conditions are assumed to be time invariant 


3.7 Unavailability of Some of the States 

Since all the states are not available, Observers are designed for reconstructing the 
unavailable states. The plant parameter matrices are to he known beforehand to design the 
Observers. But, these matrices are not known since they change with changing load 
conditions. We can use an adaptive Observer along with the state feedback adaptive 
controller. However, what happens when the two systems namely adaptive observer and 
adaptive controller designed separately are used together is not known. Another thing is 
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that the resulting system may be very complex limiting its use in practice. So, instead of 
the above app^o^iGh, we have used the parameter matrices of the reference model for state 
observation and have studied how the system behaves with such an observer incorporated. 

Since, in steady, state the plant with the feedback controller is expected to to 

tWe referev^ce rnode.l,ic \s hoped flsoJc , o-sy-mptoticcxlly 

the state observation leads to the actual state values of the plant. This observer design is 
straight forward , see for example [42}. 


3.8 Numerical Problem 
3.8.1 Using First Order Observer 

Assuming that are available and ^ is not practically available 

(which is a reasonable assumption from the practical view point), a first order observer is 
constructed using the procedure given in Appendix 3. The final equations of the whole 
system consisting of the plant, the reference model, the adaptive controller and the 
observer are given below. 

Plant: 


Reference Model: 



*2 = Vl + V3 
*3 = V3 + V4 
*4 = Vl + Vs + V4 + 


n ft 


X = X 

“l “3 


*“2 “ W ^ W 
S ” W W 


lirr-, , I 

■■Au. - .13 2212 


(3,34) 
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Controller: 


U = kjXj + + kgij + 


(3.38) 


‘'l = 


k, = T 


m 


V2 


k, = T , 7 ,x. 


T =-’ 


xn 


‘3 ^ ‘m '3^3 

^4 ~ (3.39) 

■b [ Pi(Xi - + P3(X2 - x_^) + pg(x3 - x^^) + p^(x^ - x^^) ] 

(3.40) 


Observer; 


* = * 4 * + + ‘l’'l + ‘3*2 + ‘3*4 

A 

Xg = z - (m^x^ + nijXj 

‘l = “2Vj + “'3\-'"l‘4 

‘3 = ®l + '"3'‘ie -'"2*4 
0 

*3 = “m, + ’”3\-'"3‘4 
*4 = “mg + '"3%^ + *"3%^ 


(3.41) 

(3.42) 


(3.43) 

(3.44) 


M = [mj m 3 ] = [0 0 0.0096] 

M is selected to place the eigen value of at —20 (Appendix 3). 

The system was simulated and it was found that it takes long time to settle down 
for all operating conditions. The plots of x(t) are shown for two operating conditions in 
Figs. 3.14 and 3.15. 

3.8*2 Using Second Order Observer , 

We assume that 6^ and uf^ are available^ and ’^^1' available. A second 

order observer is constructed for this situation using the procedure given in Appendix 4. 
The equations of the plant and reference model remain unchanged as in the previous case. 
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The equations of the controller and the observer are given below. 

Controller: 

U = kjXj + kjXj + kjig + 

‘‘l = '•’a Vl 

S = Vi 

''3 “ Vs 

\ V4 

= -b [ Pi(xj - + pj(x2 - x^^) 4 PjCxg - x^) 4 p^(x^ - x^^) ] 


Observer: 


(3.46) 


(3.46) 

(3.47) 


*1 “ *'62*1 “ 111^*2 ‘li*! ■*■ ‘iS^i 


*2 = ‘62*1 + V*2 + Vl + *42*2 + 


(3.48) 


Xg = Zj “ - *"2*2 

*4 =*2 "“3*1 "“4*2 

*12 ' *22 ’ *32 ' *42 ' *62 *62 defined in Appendix 4. We place the eigen values of 

(Appendix 4) at -20, -93.3772. 

Then 


' “1 "'2 ' 


■ 0 0.1178 

. "*3 ®4 . 


0 -2.6174 

s • 


(3.50) 


The systezn ia siznulated under different operating conditions and found that it does 
not pve satisfactory performance for all the operating conditions considered. The plots of 
x(t) are shown in Figs. 3.16 and 3.17. 


x3 

2 


A 

-4 


28 

1 

( 

x4 

i 

r 

( 


36 


3.8.3 Using Third Order Observer 


We assume that only is available and » e^ ^ and are unavailable. For 
this case, third order observer is constructed using the procedure given in Appendix 6. The 
equations of the plant and reference model remain unchanged as in the previous cases. The 
equations of the controller and observer are given below. 

Controller: 

u = + kgXg + + r (3.51) 

Vi 

^2 “ *^111 V2 
‘‘a = Va 

‘'4 = V4 ^ 

= -M Pl(*l - + P2(*2 - + P4(*4 - J 


(3 63) 


Observer: 


®l = "*i V’a ‘‘13“2 


*2 “ ”‘2“mj*l “2'‘inj)‘3 + + ‘23*2 

*3 = \ + \ + V® 

a 

Xj « *j - nijXj 

^ ~ *2 “ “2*2 
*4“ *3 ^“3*2 

tjg , tjj ai^ t^ we defined in Appendix B, 


(3.64) 


(3.66) 
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The system was simulated and it was found that the controller along with this third 
order observer gives very satisfactory performance over a large range of operating 
conditions and the system settles down quickly following a disturbance. The plots are 
shown in Figs. 3-lS to 3 23. 

So, from the above study we can say that the presented model reference adaptive 
PSS works very satisfactorily, when only is available, the other states being 
reconstructed using a third order observer. It should he mentioned here, that traditionally 
only the speed signed is used as input to the conventional nonadaptive PSS- 
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CONCLUSIONS 


41 Conclusions 

Fast acting, high gain modem voltage regulators often cause low frequency 
sustained oscillations in power systems. The conventional non adaptive PSS, used to damp 
out these oscillations, can not perform satisfactorily over a wide range of operating 
conditions. Modern control techniques involve use of model reference adaptive control or 
gain scheduling for the design of PSS to overcome this problem. MRAO is easier than STR 
and good adaptation speeds can be achieved. Further, this technique ensures the stability 
of the system which is lacking in STR techniques. In the present work, we have used 
MBAO technique for the design of PSS to a power system consisting of an alternator 
connected to an infinite bus through a transmission line. By simulating the overall system 
under various operating conditions, it has been shown that the PSS works satisfactorily 
over a wide remge of operating conditions, This design of PSS requires the feedback of all 
the state variables. Since, in practice, all the state variables are not avsilable, the 
unavailable states are reconstructed using observers and the performance of the system is 
studied by simulating over a wide range of operating conditions. The observers used the 
parameters of the reference model. 

It has been shown that an observer which reconstructs and using 

knowledge of only gives highly satisfactory results, over a wide range of operating 
conditions, when used in conjunction with the state feedback MRAC scheme. It should be 
noted that for conventional PSS is the commonly used input signal. Thus with the 
conventionally available instrumentation with generators, the proposed adaptation scheme 
can easily be implemented, is a significant finding of this thesis. 
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4.2 Scope for Further Work 

There is scope for extending the present work to the following areas. 

1. In practice, a power system consists of very large number of machines. So, the present 
study can be extended to multi machine systems. 

2. MRAO schemes using hyperstability approach can be developed to investigate if simpler 
and better controllers can be achieved. 

3. Here, we have used a linearised model of synchronous machine. This work can be 
extended to the design of MHAPSS using an exact nonlinear model of alternator and the 
results can be validated by simulation on a micromacbine system 

4. It can be investigated whether a MRAO scheme in conjunction with a conventional PSS 
will yield simpler and better adaptive controllers 
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APPENDIX 1 

MODEL OF POWER SYSTEM (Heffron-PMlips Model) 


The power system considered consists of a synchronous generator connected to an 
infinite bus through a transmission line as shown in Fig. A. 1.1. The well-known Eeffron 
Philips model for such a power system which represents the linear small signal system 
dynamics around an operating point can be given using standard notation, as below (see for 
example [34] ): 



Pig A. 1,1 Generator Connected to An Infinite Bus 
Through a Truismission line 



2 


Swing equations: 

= p_ A - p 




mA - *^eA 1) 

where all quantities are in p.u except 6j^ and which are expressed in radians and 
radians/sec, respectively. The coefficient is usually small and can be neglected. 
Furthermore, turbine— governor dynamics can be neglected which implies that can be 
set equal to zero. With these approximations, the power system can be represented in the 
block diagram form as in Fig. A 1.2. 

The power system Heffron Philips model can be represented in state space form as 

*a = ‘^a 

"A“( 2H 2H )®qA 


®fdA = (f) ®«A + «A + + (f) “ 

e c e e 


(Al.2) 


Defining 


“l = ( 3H ) 


- / ^ p - •) 

r *' 
e 




“4 = (^) 

-K K. 

K 

b = (^) 


(A.1.3) 


' 1.'*' 


' u' 


Then, the power system can be represented by the equation 
X = Ax + Bu 

where Jt = [ "a ®qA «hA 
0 10 0 
A« 1 ^ ° "5 ° 


^ “3 "4 

0 


a 

“6 'V J 


and 


B = [0 0 0 b]‘ 
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APPENDIX 2 

SYNCHRONOUS MACHINE AND POWER SYSTEM DATA 


The power system data considered is ^ven below (values are in PU to machine base 
whereever relevant). 

GeneratoiLEarameters : 

X^=1.14, Xq=0.66 X^ = 0.24, 

H = 1.5 sec, = 314.0 rad/sec. 

N&twor K Parameters : 

Rg = 0 024, Xg = 0 7 and = 1.0 PU. 

where 

t are the direct axis synchronous emd transient reactances respectively, 

is the quadrature axis synchronous reactance, 

is the direct axis transient open circuit time constant, 

K and r are the gain and time constants of the exciter respectively, 
e e 

H is the inertia constant of the synchronouB machine, 
o/q is the nominal synchronous speed, 

R andX are the resistance and inductive reactance of the transmission line. 
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APPENDIX 3 

CONSTRUCTION OF FIRST ORDER OBSERVER 


4 


We aasume that , e^^^ are available and is not available. Since the 

parameter matrices A and B of the plant are not known, the observer is designed using the 

parameter matrices A andB of the model. 

mm 

So 


and 



y = [xixjxj* (A.3.1) 

Since the output vector consists of Xj , and x^ , the state equations can be 


expressed in standard form, [42], as 



From [42] the observer is given by 

i e (Ajj + MAjj)^ + (Bj+MB2)u + [Ajj+MA^ — (A^+MA2j)hQy 



(A.3 2) 


(A.3,3) 

(A.3.4) 
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where M = ] is chosen feo place the eigen values of at any desired 


locations 
Here 

All = K J 
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Aj 2 - [0 0 a ], 

4 

Aoi = [0 a ]* 


”2 “7 


^22 = 


m. 


= [0] and 
Bj=[0 0b]* 


a a a 
m 

0 


6 “8 J 


(A.3.6) 


We can derive the following equations Cor the observer by substituting equation 
(A.3,5) in equations (A.3.3) and (A.3 4) 

(A.3.6) 
(A.3.7) 


^ = V + + ‘i*l + ‘2*2 + ‘3*4- 


Xg = z - (m^Xj + m^x^ +ni^^), 

t^ » nirta^ + 

1 2 2Sj 3 iQm 1 4 

1 0 

*3'=\ + ®3'‘nig““3‘4' 

‘4 = % + ® 2 \ + 


(A.3.8) 


The equations (A.3.6), (A.3.7) and (A.3.8) describe the required first order observer. 
Note that Xj is the value of the state variable e*^ of the plant aa given by the obaerver and 
Xj, Xjj and x^ are the state variables and respectively of the plant 


APPENDIX 4 

CONSTRUCTION OF SECOND ORDER OBSERVER 


6 


We assume that 6^ and are available and e^^^ and are not available. Since 
the parameter matrices A and B of the plant are not known, the observer is designed using 
the parameter matrices A and B of the model. 

So 


' ^1 ■ 
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0 ■ 


*1 ■ 


■ 0 ■ 


*2 
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“2 

0 




0 
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“4 


*3 
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0 

u 

. *4 , 


a 

L “6 

a 

“6 

a 

m? 

a 

“8 - 


. *4 


b 



and 

y = [ Xj Xj ]* (A.4.1) 

Since the output vector consists of and the state equations can be expressed 
in standard form, [42], as 


f • 
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a 
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0 ’ 
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‘ 0 ■ 
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X.. 
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4 



*1 
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“1 
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*2 
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From [42] the observer is given by 

i = (Ajj + MAj^z + (Bj+MB2)u + [Ajj+MAjj - 


(A.4.2) 


X = 



-M' 

I 


y 


(A.4.3) 

(A.4,4) 


7 


where 


M is a 2x2 matrix chosen to place the polra of at any desired locations 


^21 


M = 


Here 


Au = 




Aqi = 


A22 = 


mj 

”1 

"*3 


a 

a 

“3 

“4 

a 

a 

“7 

“8 

0 

0 

a 

a 

“6 

“^6 

0 

0 * 

a 

“2 

0 

0 

1 ■ 

a 

“1 
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®1 = [0 b]* and 

B 3 = to 0]* (A.4.6) 


Wfl can derive the following equations for the obeerver by substituting equatio 
(A.4,6) in equations (A.4.3) and (A.4.4) 


h = * 51*1 + *11*1 *11*1 

*1 = * 01*1 + V® 

*3 “ *1 ” “* 1*1 “ “‘ l *! 


(A4'fi) 
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x^=z^~ mgXj - 




Where 


*12 “ ™2“mj "*1*62 ~ 

*22 = "'l““*2*62“™4 V’ 

4 

*32 “ Vg + “ "'l*62 ” ‘”3®mg’ 

*42 ~ “mg ■'■ "*3 “ ®2*62 “ "*4“mg> 
tii 4 A = a«. + and 


o2 m,y 4 


(A.4.8) 


The equations (A.4.6), (A.4.7) and (A.4.8) describe the required second order 

A A 

observer. Here and are the state variables e^^ and e^^ as reconstructed by the 
observer. The variables x^ and x^ are the state variables and respectively of the 
plant. 



APPENDIX 5 

CONSTRUCTION OF THIRD ORDER OBSERVER 
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We assume that only 6/^ is available and e^^^ , and e^^ are not available. Since 
the pcurameter matrices A and 6 of the plant are not known, the observer is designed using 


the parameter matrices A andB of the model. 

mm 

So 


and 


■ *1 ■ 


*2 



It 
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m 3 
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(A.6.1) 


Since the output vector consists of only , the state equations can be expressed in 
standard form, [42], as 
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(A.6.2) 


From [42] the observer is given by 

i = (Ai + MAjj). + (Bi+MBj)u + (A,j+MAjj - (Aj^+MA2,)MJy 


{A.5,3) 
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where 


(A.5 4) 


M le 8> 3x1 mzitnx chosen to pl&ce the poles of ctt any desired locations. 

M = [ nij ]‘ 


Here, 




0 0 

a a 

m3 


a a 
m. m. 
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8 J 


^1 = [ “mj \ 0 1 

Ajj = [0] Bj = [0 0b]*and 

Bj=[0] (A6.6) 


We can derive the following equations for the observer by substituting the equation 
(A. 6.6) in the equations (A.5.3) and (A.6.4): 

*1 = + ‘13*2 

*2 = + ('‘m 3 + “2'‘m >2 + *23*2 

*3 = (Vj + "'3'‘m >1 + (“m^ + “3“m5>*2 + + V 2 + 

(A.6.6) 

A 

Xj = -nijXj 

A 

Xg sa Bjj - 

i^rsZj-mgXj (A-6-7) 
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where 

tj3 = 1 - ml 

‘23 - - t'"l'"2“», + "*2('‘m, + ”qO + “aVl 

i O 4 4 

‘33 = \ + “sVj) + + ”3%^ 

(A.6.8) 

The equations (A.5.6), (A.5.7) and (A.5.8) describe the required third order 

A It A 

observer. Here and are the variables and as reconstructed by the 

observer. The variable x^ is the state variable u>^ of the plant. 



